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Imaging modalities including magnetic resonance imaging and X-ray computed tomography are 
established methods in daily clinical diagnosis of human brain. Clinical equipment does not provide 
sufficient spatial resolution to obtain morphological information on the cellular level, essential for applying 
minimally or non-invasive surgical interventions. Therefore, generic data with lateral sub-micrometer 
resolution have been generated from histological slices post mortem. Sub-cellular spatial resolution, lost in 
the third dimension as a result of sectioning, is obtained using magnetic resonance microscopy and micro 
computed tomography. We demonstrate that for human cerebellum grating-based X-ray phase tomography 
shows complementary contrast to magnetic resonance microscopy and histology. In this study, the 
contrast-to-noise values of magnetic resonance microscopy and phase tomography were comparable 
whereas the spatial resolution in phase tomography is an order of magnitude better. The registered data with 
their complementary information permit the distinct segmentation of tissues within the human cerebellum. 

Minimally or non-invasive surgical interventions on the brain are generally based on less detailed patient- 
specific imaging (pre- and intra-operative) and high-resolution generic anatomical data. The generic 
knowledge of brain microstructure currently relies on two-dimensional (2D) evaluation of histological 
slices, which requires time-consuming sectioning, staining, and (manual, expert-based) image analysis. The 
advantages of histology are sub -micrometer resolution and excellent tunable contrast depending on the staining 
protocol 1 ' 2 . However, the technique does not provide isotropic three-dimensional (3D) information of the soft 
and hard tissues, as is regularly obtained using magnetic resonance imaging (MRI) and X-ray computed tomo- 
graphy (CT). Furthermore, it is hard to imagine image-guided neurosurgery 3 without MR guidance. 
Conventional MRI is a well-established technique for brain imaging, and yields superb contrast between white 
and grey matter but only limited spatial resolution. At present, medical MR systems produce images with typically 
sub -millimetre voxel lengths. Small animal MR microscopy (uMRI) scanners equipped with stronger gradient 
systems and operating at higher magnetic field strengths have been used to visualize mouse brains in vivo 4,5 and 
pieces of human brain post-mortem 6,7 with voxel sizes of a few tens of micrometers. Conventional CT also 
provides fully quantitative 3D data ultimately reaching higher spatial resolution than MRI. The CT contrast 
for brain tissue, however, is weak. Despite low inherent brain tissue contrast, micro-CT (uCT) as a non-destruct- 
ive technique has already been used to correct deformations induced in data sets obtained by histological 
sectioning during the preparation of histology 8 . In order to obtain a better contrast than in these absorption- 
based studies, synchrotron radiation-based phase-contrast uCT (PC-uCT) 9,10 was recently applied to human 
brain 11 . X-ray phase-contrast methods are based on the phase shifts of X-ray waves penetrating the specimen 
which is related to the decrement 5(x,y,z) of the real part of the refractive index distribution. A variety of PC-uCT 
methods based on crystal interferometry 12 , propagation-based (or inline) contrast 1314 or analyzer-based imaging 
(ABI) 1516 are available today. More recently, grating interferometry 17 " 19 (also known as Talbot interferometry) 
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Figure 1 | Cerebellum block extraction. The photograph (right) and MR-slice acquired using a Verio 3T whole body scanner (Siemens Health Care, 
Erlangen, Germany) illustrate the location and size of the cerebellum within the donated body. 



was developed, which reaches the especially high contrast required 
for identifying structures with small differences in electron density as 
is the case in brain 11 ' 20 ' 21 . 

In the light of these inherent advantages and disadvantages of 
MR- and CT-based techniques and histology, the aim of the present 
communication is to elucidate how far the combination of the com- 
plementary methods |iMRI, grating-based PC-|aCT and histology 
yields additional information on the microanatomy of the human 
cerebellum. 

Results 

3D imaging. A specimen taken from the human cerebellum (Figure 1) 
was visualised in 3D using grating-based PC-|iCT and |iMRI. Figure 2 
shows virtual cuts through identical regions of the brain tissue obtained 
with the two imaging techniques after non-rigid 3D registration 22 . 
In the mixed T1/T2* -weighted |iMRI images white matter structures 
appear darker than grey matter structures. Grey matter consists of two 
components: the stratum granulosum and the stratum moleculare. 
Because of the weak MR-contrast between these layers, however, it is 
impossible to segment these structures with an intensity-based algo- 
rithm. PC-|iCT, on the contrary, provides high contrast between these 
grey matter layers: with the stratum granulosum appearing as a bright 
feature and the stratum moleculare as the dark outer structure. The 
white matter (interior dark region in PC-|iCT) exhibits ^-values very 
similar to the stratum moleculare and can only be segmented because 
the stratum granulosum is located between the two tissues. Therefore, 
an intensity-based segmentation is also impossible in this case. This can 
be deduced from the related histograms of the 3D data in Figure 2. The 
histogram of the |iMRI data (green-coloured triangles in the figure) 
shows a shoulder (left) associated with white matter, and a peak asso- 
ciated with grey matter (right). The histogram of the PC-|iCT data (red- 
coloured squares) contains two peaks and a shoulder. The first peak 
corresponds to the stratum moleculare, which was in direct contact with 
formalin during the whole fixation period of the brain. The second peak 
is a superposition of the stratum moleculare that was in contact with 
formalin only for restricted periods of time and the white matter. The 
shoulder on the right originates from the stratum granulosum. 



Tissue contrast. The images in Figure 2 acquired with PC-|iCT and 
|iMRI show the characteristic morphology of the human cerebellum. 
The contrast-to-noise ratio c between two anatomical structures can 
be defined as 



\Xc\-Xc2\ 

FWHM back 



(1) 



where x ci are the centre positions of the related peaks in the histo- 
grams and FWHM back is the full-width-at-half-maximum of the 
background signal, here formaldehyde solution and perfluoropo- 
lyether, respectively. This definition of the denominator was chosen 
for practical reasons (see supplementary material). Using Equation 1, 
six c-values between the four anatomical structures were determined 
for PC-|iCT and jiMRI (Figure 3). The highest contrast in the PC- 
|iCT images is the contrast between stratum granulosum and tissue 
in contact with formalin (E) and has the value of c PC -^cT_max = 6.4 ± 
0.2. Using |iMRI a maximum value of c flMRI max = 9.5 ±0.1 was 
determined, corresponding to the contrast between white matter and 
tissue in contact with formalin (C). 

Joint histogram. The superposition of peaks in the individual 
histograms prevents the segmentation of the anatomical features. 
In particular, it is impossible to differentiate between stratum 
moleculare and stratum granulosum in the |iMRI data, and 
between white matter and stratum moleculare in the PC-|iCT data. 
Combining |iMRI and PC-|iCT using a joint histogram, however, 
allows the distinct discrimination between the characteristic soft 
tissues (Figure 4). To do so, the first step was the identification of 
the 2D maxima in the joint histogram corresponding to the tissues of 
interest using a multi-Lorentzian fit. Subsequently, the intersections 
of the Lorentzians were calculated. These intersections characte- 
rize the interfaces between the anatomical structures and allow a 
distinct differentiation between white matter, given in red, stratum 
granulosum in blue, and stratum moleculare in yellow. In addition, 
the green- coloured peak indicates regions of stratum moleculare, 
which were in direct contact with formalin during the whole 
fixation period. 
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Figure 2 | X-ray phase micro computed tomography and MR microscopy of the cerebellum. The selected orthogonal virtual cuts through the specimen 
and the corresponding histograms of the 3D data illustrate the power and complementary nature of the two techniques for imaging the human 
cerebellum. 



Spatial resolution. In order to evaluate the spatial resolution of 
|iMRI, PC-|iCT, and histology, firstly line plots through the 
registered data of the human cerebellum (see Figure 5) were 
inspected. The spatial resolution of the imaging techniques can be 
deduced by comparing the roughness of the line plots. Histology has 
the best spatial resolution followed by PC-|iCT then |iMRI. To 
quantify this observation, the spatial resolution X S r was calculated 
by Fourier analysis 11 ' 23 . The radial spectral power (rSP) of a structure 
of interest in a 100 X 100 pixel region of a selected slice (rSP struc ) was 
calculated and divided by the rSP of a background region of same size 
(rSP back ) (i.e. of formalin for PC-|iCT, perfluoropoly ether for |iMRI 
and glass without brain tissue for histology). This ratio was plotted 




A B C D E F 

Tissue relations 

Figure 3 | Tissue contrast in human cerebellum. The above data compare 
contrast values between uMRI and PC-uCT for white matter and stratum 
granulosum (A), white matter and stratum moleculare (B), white matter 
and tissue in contact with formalin (C), stratum granulosum and stratum 
moleculare (D), stratum granulosum and tissue in contact with formalin 
(E) and between stratum moleculare and tissue with formalin (F). 



against the spatial frequency. The peak-to-peak distance of the first 
frequency larger than twice the mean value of the baseline (noise) 
was divided by two to determine the spatial resolution 24 . The X SR - 
values corresponded to (6 ± 1) |im for histology, to (20 ± 5) |im for 
PC-uGT, and to (210 ± 40) Jim for uMRI. 

Multimodality. The morphology of the cerebellum can be properly 
visualized combining (xMRI, PC-|iCT, and histology. Figure 6 
illustrates one possible way of combining the 2D images of the 
three techniques used in this study. Each technique is represented 
by one of the colour channels red/green/blue (RGB) (see Figure 6, left 
image). Such a colour representation directly illustrates the contri- 
butions of the individual techniques. For example the yellow colour, 
which indicates the stratum granulosum, is mainly the superposition 
of contrast from the PC-(iCT and |iMRI, whereas the purple colour 
(white matter) primarily arises from the PC-|iCT and the histology 
signals. 

Discussion 

Grating interferometry is based on the detection of the deflection 
angles of the X-rays passing through the specimen. After reconstruc- 
tion, a 3D distribution of the decrement 5 of the real part of the 
refractive index of the specimen can be obtained, which is directly 
related to the electron density distribution 



r X 

d(x,y,z) = ^p e (x,y,z) 



(2) 



with the classical electron radius r e and the X-ray wavelength X 25 . 

In MRI, contrast between different tissues is determined by their 
relative spin densities p s , their characteristic relaxation properties at 
the purported magnetic field strength, i.e. their distinguished 7\ 
(longitudinal) and T 2 resp. T 2 * (transversal) relaxation times, and 
the imaging sequence and parameters used (in particular the echo 
time TE, the repetition time TR and the flip angle a). In a very general 
manner, the measured amplitude of the signal for a gradient echo 
sequence can be described as 
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Figure 4 | Segmentation by means of a joint histogram. From the joint histogram of the phase tomography and MR microscopy data it is possible to 
distinctly segment the stratum granulosum, stratum moleculare, white matter, and brain tissue that was in direct contact with formalin for a longer period 
of time than other parts. The white arrows indicate the positions of the related line plots. The colours of the virtual cuts of PC-uCT and uMRI on the right 
side correspond to the colours assigned to the four peaks in the joint histogram. 
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Figure 5 | Spatial resolution of histology, X-ray phase tomography, and MR microscopy. Comparable line plots through selected slices demonstrate that 
the spatial resolution decreases from histology via PC-uCT to uMRI. 
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Figure 6 | Multimodal imaging of the cerebellum. The fusion of slices from PC-uCT, uMRI, and histology using the RGB channels directly visualizes the 
contributions from the three imaging techniques. The combination of the three techniques permits the distinction of otherwise indistinguishable 
anatomical features, such as the separation between white matter (violet), stratum granulosum (yellow), stratum moleculare (green) and blood vessels of 
different sizes (purple and white arrows). The enlarged regions (with different magnification) show the light blue coloured small capillaries. 



S(p,Ti,T 2 ,TE,TR) oc p s sin 



oc(l-e — 



where 



l - = JL +yA£ 0 with yAB 0 ^60Hz 



(3) 



(4) 



takes into account the refocusing of magnetization components that 
were dephased under the readout gradient of the imaging sequence, 
while dephasing of magnetization resulting from static field inho- 
mogeneities such as susceptibility effects contribute to the image 
contrast. In order to optimize tissue contrast, dedicated imaging 
sequence parameters were chosen based on literature values for grey 
and white matter at magnetic fields of 9.4 T 26 . 

Before the acquired data can be compared and combined, they 
must be registered with each other in 3D space. It was shown that the 
3D registration using the maximization of the mutual information 22 
gives precise results for the registration of tomography data from 
different imaging modalities 27 . A rigid registration including three 
degrees of translation and three degrees of rotation seems to be 
sufficient for the 3D data set acquired by PC-|iCT and |iMRI. 
Although the brain tissue undergoes shrinkage during formalin fixa- 
tion, further deformations over long-time periods are minimal. 
Nevertheless, an affine registration with an additional three scaling 
factors compensates for the error in individual pixel sizes. In our case 
a non-rigid registration was required to correlate the corresponding 
voxels between the imaging modalities. This approach accounts for 
possible deformations within the brain tissue induced during the 
transfer of the specimen from the cylindrical polyethylene container 
(filled with formaldehyde) used in PC-|iCT measurements into the 
container (with perfluoropoly ether) for jjMRI. Because of the defor- 
mations induced during sectioning and staining, the histological slice 
was also non-rigidly registered with the tomography data. 

The virtual cuts through the |iMRI and PC-|iCT data sets in 
Figure 2 elucidate the pros and cons of each method. It is well known 
that MRI provides superb contrast between white and grey matter. 
White matter (dark region) can easily be distinguished from the grey 



matter (bright region) in our mixed T1/T2* -weighted |iMRI images. 
On the contrary, differentiation between layers of grey matter tissues, 
in particular between stratum granulosum and stratum moleculare, 
is hardly possible. Here, PC-|iCT provides a clear distinction between 
these anatomical structures. White matter (inner dark region), how- 
ever, has a comparable electron density to stratum moleculare (outer 
dark region). These observations are quantified using Equation 1 and 
summarized in Figure 3. With the definition of the contrast-to-noise 
ratio c presented in this manuscript, a value c > 1 is required to 
distinguish anatomical structures using intensity-based algorithms. 
Therefore, one can conclude that the PC-(iCT contrast is sufficiently 
high to distinguish between the four tissues whereas |iMRI contrast is 
restricted to white and grey matter. Additionally, the PC-|iCT pixel 
size is a factor of nine smaller than jaMRI, although the maximal 
c- value observed in the PC-|iCT results is slightly lower than the 
maximal |iMRI contrast-to-noise ratio. Binning the tomography 
data reduces the spatial resolution but improves the contrast 28 . 
Therefore, we can state that the contrast-to-noise between the two 
methods yields comparable values. 

In addition to the contrast-to-noise ratio c, we also consider the 
spatial resolution X SR to assess image quality. Following the work 
of Thurner et al. 28 , we combine these two values in the quality factor 



c 

^SR 



(5) 



Using the maximal contrast values of PC-|iCT and |iMRI and the 
corresponding spatial resolution values, one obtains quality factors 
of (0.32 ± 0.09) urn" 1 for PC-uCT and (0.05 ±0.01) urn" 1 for uMRI. 
Based on this definition, PC-|iCT yields six times better images than 
|iMRI. More importantly, however, the complementary results of the 
two data sets can be derived using a joint histogram as presented in 
Figure 4. Three anatomical structures can be segmented, i.e. white 
matter (red), stratum granulosum (blue) and stratum moleculare 
(yellow). Furthermore, due to its high sensitivity, PC-|iCT allows 
the detection of differences within the tissue induced by different 
periods of formalin fixation. The green regions in Figure 4 cor- 
respond to domains of stratum moleculare, which were in direct 
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contact with formalin for the whole fixation period of 120 days and 
not only for 20 days. 

In summary, we demonstrate the advantages and disadvantages of 
PC-|iCT, uMRI and histology and present a method to combine data 
from these complementary techniques in order to segment the ana- 
tomical layers of the human cerebellum post mortem. The clear 
segmentation of the three layers allows us to quantify the ratios 
between these volumes. Given that cell death in the cerebellum is a 
common cause of neurological disease, e.g. in cerebellar ataxias 29 and 
that abnormalities in the cerebellum have more recently also been 
implicated in autism and neuropsychiatric disorders 30 ' 31 methods for 
a reliable quantitative assessment of cerebellar morphology and 
pathology are certainly needed. Unfortunately, the currently avail- 
able stereological methods for the cerebellum are very time consum- 
ing and often lack precision 31,32 . The presented method has the 
potential to significantly reduce the effort and improve the quality 
of quantitative stereological assessments of the cerebellum. 

Methods 

Specimen preparation. The brain from the donated body of a 68-year-old male was 
extracted at the Institute of Anatomy (University of Basel, Switzerland) within 48 h 
after death. The whole brain was transferred to 10% formalin solution for fixation. All 
procedures were conducted in accordance with the Declaration of Helsinki and 
according to the ethical guidelines of the Canton of Basel. The small cerebellum block 
(6 X 6 X 11 mm 3 ) was extracted at the University Hospital Zurich, Switzerland, after 
approximately three months of fixation (Figure 1). During the measurements the 
specimen was placed in a 0.5 mL Eppendorf container filled with 4% formalin 
solution for PC-uCT and with perfluoropolyether (Fomblin® Y-LC 80, Solvay Solexis, 
Bollate, Italy) for uMRI in order to reduce the background signal and susceptibility 
artefacts at the tissue-liquid interface 33 . 

Grating-based X-ray phase tomography. A detailed description of grating 
interferometry can be found in the literature 1119 ' 34 ' 35 . The experiment was performed 
at the beamline ID 19 (ESRF, Grenoble, France) at the interferometer installed there 36 
with X rays from the U32 undulator with a gap set to 15.15 mm. Using a double- 
crystal Si(lll)-monochromator in Bragg geometry a photon energy of 23 keV was 
selected. The imaging of the entire specimen required two scans to obtain a field of 
view 11 mm high and 6 mm wide. The beam-splitter grating^ (made at the Paul 
Scherrer Institut, Villigen, Switzerland) with a periodicity of pi = 4.785 um had a Si 
structure height of 29 um, which corresponds, within a few percent, to a phase shift of 
7i for 23-keV X rays. The analyzer grating g 2 (made at the Karlsruhe Institute of 
Technology, Eggenstein-Leopoldshafen, Germany) had lines of gold with a 
periodicity ofp 2 = 2.400 um and 50 um structure height. This corresponds to an 
intensity transmission of only about 6 X 10" 3 through the teeth of the analyzer grating 
comb, i.e., a contrast of 100%. With a distance between source and interferometer of 
150 m and a distance between the gratings of d = 479.4 mm (corresponding to the 9 th 
Talbot order), the ratio of the grating periods p 2 /pi was matched to the beam 
divergence 35 . The sample container was immersed in a water tank with parallel 
polymethylmethacrylate plates in order to minimize artefacts due to X-ray phase 
curvature induced by a conical container-air interface. The specimen was located 
10 cm upstream of the beam-splitter grating. The detector was a lens- coupled 
scintillator and charge-coupled device (CCD) system using a FreLoN 2K (Fast- 
Readout, Low-Noise, ESRF Grenoble, France) CCD with 2048 X 2048 pixels and an 
effective pixel size of 5. 1 um, placed approximately 3 cm downstream of the analyzer 
grating. Projection radiographs were taken in 1501 steps over a range of 360°. At each 
projection angle, four phase- stepping images with an exposure time of 1 s each were 
taken over one period of the interference pattern. 

Magnetic resonance microscopy. The uMRI measurements were carried out on a 
9.4 T, 30 cm horizontal small animal MR unit (Bruker BioSpec, Bruker BioSpin MRI, 
Ettlingen, Germany) equipped with a transceiver cryogenic quadrature radio 
frequency (RF) surface coil (CryoProbe™) designed for MRI of the mouse brain at 
increased sensitivity. High-resolution structural images were acquired with a T 2 *- 
weighted 3D FLASH sequence with isotropic voxel size of 45 um (field of view (FOV) 
1.35 X 0.90 X 1.35 mm 3 , acquisition matrix (MTX) 300 X 200 X 300, echo time (TE) 
12 ms repetition time (TR) 400 ms, flip angle (a) 15°, number of averages (NA) 1 at a 
total acquisition time of 400 minutes. 

Histology. Sections of 50 um thickness of the cerebellum block were prepared using a 
microtome with vibrating blade HM 650 V (Microm International GmbH, Walldorf, 
Germany). For the staining procedure, the sections were mounted on gelatinized 
slides and stained for Nissl with cresyl violet 1 . Microphotographs of selected slices 
were taken using a Leica MZ16 microscope and DFC420-C digital camera. 

Data processing. The reconstruction of the phase-contrast dataset was carried out 
using a modified filter kernel (Hilbert transform) in combination with a standard 
back-projection algorithm 34 ' 37 ' 38 . Because the signal intensity in the uMRI images 



varies with the coil sensitivity profile, further processing was performed. For the 
correction of the inhomogeneities a modified fuzzy C-means algorithm was applied 39 . 
For a reasonable analysis of the joint histogram, the 3D datasets were non-rigidly 
registered using the classical maximization of mutual information (MI) principle 40 ' 41 . 
The chosen registration algorithm is based on the adaptive hierarchical image 
subdivision strategy, which decomposes the non-rigid matching problem into 
numerous local affine registrations of image segments of decreasing size 22 . The local 
registrations were again performed using the classical maximization of MI principle. 
The local registration parameters were identified using the Powell multi-dimensional 
search algorithm 42 such that the MI between the reference and the floating sub -images 
was maximized. The hierarchical image splitting strategy 43 recently underwent 
several improvements and was extended to 3D 22,44 . For region selection in the joint 
histogram, the line plots were approximated with a multi-Lorentzians fit using the 
Levenberg-Marquardt algorithm in OriginPro 7.5 (OriginLab Corporation, 
Northampton, USA). 



1 . Morel, A. Stereotactic Atlas of the Human Thalamus and Basal Ganglia. (Informa 
Healthcare New York, 2007). 

2. Schaltenbrand, G. & Wahren, W. Atlas for Stereotaxy of the Human Brain. 
(Thieme, Stuttgart, 1977). 

3. Jolesz, F. A., Nabavi, A. & Kikinis, R. Integration of interventional MRI with 
computer-assisted surgery. /. Magn. Reson. Imaging 13, 69-77 (2001). 

4. Benveniste, H. & Blackband, S. MR microscopy and high resolution small animal 
MRI: Applications in neuroscience research. Prog. Neurobiol. 67, 393-420 (2002). 

5. Ahrens, E. T. et al. MR microscopy of transgenic mice that spontaneously acquire 
experimental allergic encephalomyelitis. Magn. Reson. Med. 40, 119-132 (1998). 

6. Fatterpekar, G. M. et al. Cytoarchitecture of the human cerebral cortex: MR 
microscopy of excised specimens at 9.4 Tesla. Am. J. Neuroradiol. 23, 1313-1321 
(2002). 

7. Fatterpekar, G. M. et al. MR microscopy of normal human brain. Magn. Reson. 
Imaging Clin. N. Am. 11, 641-653 (2003). 

8. Germann, M. et al. Strain fields in histological slices of brain tissue determined by 
synchrotron radiation-based micro computed tomography. /. Neurosci. Methods 
170 (1), 149-155 (2008). 

9. Fitzgerald, R. Phase- sensitive X-ray imaging. Rhys. Today 53 (7), 23-26 (2000). 

10. Momose, A., Fujii, A., Kadowaki, H. & Jinnai, H. Three-dimensional observation 
of polymer blend by X-ray phase tomography. Macromolecules 38 (16), 7197- 
7200 (2005). 

11. Schulz, G. et al. High-resolution tomographic imaging of a human cerebellum: 
Comparison of absorption and grating-based phase contrast. /. JR. Soc. Interface 7, 
1665-1676 (2010). 

12. Bonse, U. & Hart, M. An X-ray interferometer. Appl. Rhys. Lett. 6 (8), 155-156 
(1965). 

13. Snigirev, A. A., Snigireva, I., Kohn, V., Kuznetsov, S. & Schelokov, I. On the 
possibilities of x-ray phase contrast microimaging by coherent high-energy 
synchrotron radiation. Rev. Sci. Instrum. 66 (12), 5486-5492 (1995). 

14. Cloetens, P., Barrett, R., Baruchel, J., Guigay, J.-P. & Schlenker, M. Phase objects in 
synchrotron radiation hard x-ray imaging. /. Rhys. D: Appl. Rhys. 29, 133-146 
(1996). 

15. Forster, E., Goetz, K. & Zaumseil, P. Double crystal diffractometry for the 
characterization of targets for laser fusion experiments. Krist. Tech. 15 (8), 
937-945 (1980). 

16. Davis, T. J., Gao, D., Gureyev, T. E., Stevenson, A. W. & Wilkins, S. W. Phase- 
contrast imaging of weakly absorbing materials using hard X-rays. Nature 373, 
595-598 (1995). 

17. David, C, Nohammer, B., Solak, H. H. & Ziegler, E. Differential X-ray phase 
contrast imaging using a shearing interferometer. Appl. Rhys. Lett. 81, 3287-3289 
(2002). 

18. Momose, A. et al. Demonstration of X-ray Talbot interferometry. Jpn. J. Appl. 
Rhys., Part 2 42 (7B), 866-868 (2003). 

19. Weitkamp, T. et al. X-ray phase imaging with a grating interferometer. Opt. 
Express 13, 6296-6304 (2005). 

20. Pfeiffer, F. et al. High- sensitivity phase-contrast tomography of rat brain in 
phosphate buffered saline. /. Rhys. Conf. Ser. 186, 012046 (2009). 

21. Pfeiffer, F. et al. High-resolution brain tumor visualization using three- 
dimensional X-ray phase contrast tomography. Rhys. Med. Biol. 52 (23), 
6923-6930 (2007). 

22. Andronache, A., von Siebenthal, M., Szekely, G. & Cattin, P. Non-rigid 
registration of multi-modal images using both mutual information and 
cross-correlation. Med. Image Anal. 12 (1), 3-15 (2008). 

23. Modregger, P., Liibbert, D., Schafer, P. & Kohler, R. Spatial resolution in Bragg- 
magnified X-ray images as determined by Fourier analysis. Rhys. Status Solidi A 
204, 2746-2752 (2007). 

24. Wang, L., Ho, P. P., Liu, C, Zhang, G. & Alfano, R. R. Ballistic 2-D imaging 
through scattering walls using an ultrafast optical Kerr gate. Science 253, 769-771 
(1991). 

25. Als-Nielsen, J. & McMorrow, D. Elements of Modern X-ray Physics. (John Wiley & 
Sons, Chichester, 2011). 

26. Kuo, Y.-T., Herlihy, A. H., So, P.-W., Bhakoo, K. K. & Bell, J. D. In vivo 
measurements of Tl relaxation times in mouse brain associated with different 



SCIENTIFIC REPORTS | 2 : 826 | DOI: 1 0.1 038/srep00826 



6 



modes of systemic administration of manganese chloride. /. Magn. Reson. 
Imaging 21, 234-239 (2005). 

27. Miiller, B. et al. Three-dimensional registration of tomography data for 
quantification in biomaterials science. Int. J. Mater. Res. 103, 242-249 (2012). 

28. Thurner, P., Beckmann, F. & Miiller, B. An optimization procedure for spatial and 
density resolution in hard X-ray micro -computed tomography. Nucl. Instrum. 
Methods Rhys. Res. B 225, 599-603 (2004). 

29. van Gaalen, J. & van de Warrenburg. B.P.C. A practical approach to late-onset 
cerebellar ataxia: putting the disorder with lack of order into order. Pract. Neurol. 
12, 14-24 (2012). 

30. Villanueva, R. The cerebellum and neuropsychiatric disorders. Psychiatry Res. in 
press, published online before print 22 March (2012). 

31. Andersen, K., Andersen, B. B. & Pakkenberg, B. Stereological quantification of the 
cerebellum in patients with Alzheimer's disease. Neurobiol. Aging 33, 197. 
elll-197.el20 (2010). 

32. Agashiwala, R. M., Louis, E. D., Hof, P. R. & Perl, D. P. A novel approach to non- 
biased systematic random sampling: A stereologic estimate of Purkinje cells in the 
human cerebellum. Brain Res. 1236, 73-78 (2008). 

33. Smith, B. R., Huff, D. S. & Johnson, G. A. Magnetic resonance imaging of embryos: 
An internet resource for the study of embryonic development. Comp. Med. Imag. 
Graph. 23, 33-40 (1999). 

34. Pfeiffer, F., Kottler, C, Bunk, O. & David, C. Hard X-ray phase tomography with 
low-brilliance sources. Phys. Rev. Lett. 98, 108105 (2007). 

35. Weitkamp, T., David, C, Kottler, C, Bunk, O. & Pfeiffer, F. Tomography with 
grating interferometers at low-brilliance sources. Proc. SPIE 6318, 63180S (2006). 

36. Weitkamp, T. et al. Recent developments in X-ray Talbot interferometry at ESRF- 
ID19. Proc. SPIE 7804, 780406 (2010). 

37. Faris, G. W. & Byer, R. L. Three-dimensional beam-deflection optical tomography 
of a supersonic jet. Appl. Opt. 27, 5202-5212 (1988). 

38. Pfeiffer, F., Bunk, O., Kottler, C. & David, C. Tomographic reconstruction of 
three-dimensional objects from hard X-ray differential phase contrast projection 
images. Nucl. Instrum. Methods Phys. Res. Sect. A 580, 925-928 (2007). 

39. Ahmed, M. N., Yamany, S. M., Mohamed, N., Farag, A. A. & Moriarty, T. A 
modified fuzzy C-means algorithm for bias field estimation and segmentation of 
MRI data. IEEE Trans. Med. Imaging 21, 193-199 (2002). 

40. Maes, F., Collignon, A., Vandermeulen, D., Marchal, G. & Suetens, P. Multi- 
modality image registration by maximization of mutual information. 
Mathematical Methods in Biomedical Image Analysis. IEEE, 14-22 (1996). 

41. Viola, P. & Wells, W. M. Alignment by maximization of mutual information. 
Proceedings of the fifth international conference on computer vision, 16-23 (1995). 

42. Press, W. H., Flannery, B. P., Teukolsky, S. A. & Vetterling, W. T. Numerical 
Recipes in C - The Art of Scientific Computing. (Cambridge University Press, 
1988). 

43. Likar, B. & Pernus, F. A hierarchical approach to elastic registration based on 
mutual information. Image Vis. Comput. 19 (1-2), 33-44 (2001). 



44. Schulz, G. et al. Three-dimensional strain fields in human brain resulting from 
formalin fixation. /. Neurosci. Methods 202, 17-27 (2011). 



Acknowledgements 

The authors gratefully acknowledge M. Miiller- Gerbl and P. Zimmermann for the 
organization and carrying out the extraction of the human brain, J. Kapfhammer for 
anatomical and medical information, A. Morel for the brain specimen blocking, M. Imholz 
for the histology of the specimen, S. Rutishauser for the preparation of the beam-splitter 
grating, M. Rudin for giving us the opportunity to perform the magnetic resonance 
microscopy experiments and A. Andronache for making the non-rigid registration tool 
available. The authors thank M. N. Holme for improving the readability and for correcting 
the language. T. W. acknowledges support from the French research networks (reseau 
thematique de recherche avancee, RTRA) "Digiteo" and "Triangle de la Physique" (grants 
2009-034T and 2009-79D). F. P. acknowledges financial support through the DFG Cluster 
of Excellence Munich-Center for Advanced Photonics and the European Research Council 
(FP7, Starting Grant No. 240142). The project was partially funded by the Swiss National 
Science Foundation (CR23I2_125 406) and was supported by the ESRF (proposal MD-407) 
through beam time allocation. This work was partly conducted with the support of the 
Karlsruhe Nano Micro Facility (KNMF), a Helmholtz Research Infrastructure at Karlsruhe 
Institute of Technology (KIT). 

Author contributions: 

G.S. prepared the specimen for the measurements; F.P., T.W. and CD. developed the 
phase-contrast set up; G.S., I.Z., T.W. and B.M. performed phase-contrast tomography 
experiments; C.W. performed magnetic resonance microscopy experiments; G.S. and B.M. 
analyzed the data and prepared the figures; G.S. and B.M. wrote the initial manuscript; all 
authors discussed the results and contributed to the final manuscript. 

Additional information 

Supplementary information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 
License: This work is licensed under a Creative Commons 

Attribution-NonCommercial-NoDerivs 3.0 Unported License. To view a copy of this 
license, visit http://creativecommons.Org/licenses/by-nc-nd/3.0/ 

How to cite this article: Schulz, G. et al. Multimodal imaging of human cerebellum - 
merging X-ray phase microtomography, magnetic resonance microscopy and histology. 
Sci. Rep. 2, 826; DOI:10.1038/srep00826 (2012). 



SCIENTIFICREPORTS | 2 : 826 | DOI: 1 0.1 038/srep00826 



7 



